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Andrzej DĄBROWSKI 


Effects of Surface 'Hetetogeneity in Adsorption from Non-ideal Binary Liqui 
Mixtures on Solids 


Wplyw heterogeniezńości adsorbentu na proces adsorpcji z nieidealnych dwuskładnikowych | 
, —" nieelektrolitów na powierzchniach ciał stałych - a 


BunuaHne rereporenHocTW a1COpÓeHTa Ha npouecc arcop6usi M3 ÓnHapHbIX HeM1EA.IbE 
DacTBOpOB He5JIeKTpOJIMTOB Ha NOBEPXHOCTAX THEPALIX ren 1 


INTRODUCTION 


The influence of surface 5 on the mechaniem | 
of the adsorption process at solid/liquid interface haa deen ji. 
intensively studied during last decade [ 1-15 J. Recently, a 
comprehansive survey of various approaches to liquid adsorp- bs 
tion onto heterogeneous adsorbents nas been given by Da — 
browski et al. [14]. In the cited papers both mono- 
layer as well as multilayer models of adsorbed pnase form- 
ing on heterogeneous solids were discussed. however, the it 
monograph C 14] pointed out inadequate investigations. deal- * 
ing with the role of surface heterogeneity of the adsorbent — : 


3 w = O eaa ̃ — . . ... ee ̃ DN . — 
~g - v = ` * i 
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when a liquid mixture is non-ideal in both bulk and surface 
shaban. The recent papers of Dabrowski et al. [3] 
and Da browski [8-41] do not explain this problem 
satisfactorily. On the other hand, Everett [15] 


UxvettJenoyer [16] and Hansen [17] suggested 
that the heterogeneity of the surface can play a dominant role 


as a source of imperfection in the adsorbed region. Simultaneously, 


ought to expect that (except for the adsorption onto extramely 
strongly heterogeneous surfaces where a dominant role of the 
factor in question is beyond doubt proved [5, 11, 18 ]) an 
influence of the non-ideality of the bulk and (or) surface 
phases on the mechanism of adsorption at solid/liquid interface 


is like to that arising from surface heterogeneity of the adsor- 


bent used. Hence, from theoretical view point a double possibility 


to explain of deviations of experimental results. from simplest 


form of Everett's equation (this equation reflects so-called 


ideal behaviour of adsorption system ) appears: 1. either in 
terms of so-called surface activity coefficients, or 2. in terms 


of surface heterogeneity effects. Fortunately, bulk activity 


coefficients are independent on the adsorption experiment and, 


in view of it, a role of interactions at free solution in liquid 


_ adsorption can be easily established. The determination 


separately of an influence of intermolecular interactions among 


admolecules and surface heterogeneity of the solid on imperfec- 


tion of the adsorbed phase is more complicated. 


To look for a solution of the above problem, so-called 
Langmuir-Freundlich (LF) equation will be handled by us. Pre- 
viously this type of isotherm has been introduced in liquid ad- 


| 
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sorption by means of kinetic approach by Jaroniec [19], 
however, Dabrowski etal. [27] explainea it on the ba- 
sis of both the Stieltjes transform and condensation approxima- 
tion method [20] It has been shown [3 ] that (LF) equation 
describes the monolayer adsorption from solutions consisting of 
molecules of equal sizes and this one corresponds to the symmetri- 
cal quasi~gaussian energy distribution. 

| In this paper on the basis of (LF) equation and assuming 
the simple statistical model for regular solutions, analytical 
expressions for isosteric heats of adsorption are derived. Next, 
the combined effects of surface heterogeneity and non-ideality 
of both bulk and surface phases as well as interactions among 
molecules belonging to those on the excess adsorption isotherms 
and isosteric heats of adsorption, are studied numerically. 
Besides, criteria for the formation of S-shaped excess isotherms 
are established. Next, the problem of surface activity coeffi- 
cients is widely discussed. In terms of (IF) isotherm, the 
suitable equations are proposed with the aim to separate an 
influence of heterogeneity effects and intermolecular inter- 
actions between admolecules on the imperfection of the adsorbed 
phase. Above-mentioned equations are studied theoretically and 
some literature experimental data are applied to their verifi- 
cation. It is proved that experimental results explained 
previously in terms of negative deviations of adsorbed solution 
from Raoult’s law can be explained in terms of a more realistic 
assumption respecting surface heterogeneity of the solid., 
Moreover, by means of suitable experimental data a temperature 
dependence of the heterogeneity parameter appearing in - (123 


. equation is discussed too. 
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Regardless of the fact that the crude assumptions underly 

_ (LF) equation(i.e., monolayer model of surface phase and symme- 
trical quasi-gaussian energy distribution corresponding to it) 
expression in question seems to be one of the most popular and 
useful equations used for characterizing many simple adsorption 
Systems [5-5, 21-24 |. Hence, we believe that our considerations 
| Such as below can provide valuable informations about the nature 
.of the adsorption process at solid/liquid interface. 


GENERAL CONSIDERATIONS 


= 


one of the most fundamental equations dealing with the 
liquid adsorption on heterogeneous solid surfaces, characterized 
by a continuous energy distribution has the following form C14]: 


x (7) s d 4, 104, €12) A 643986 42 (1) 
4 d 
— „here X t denotes the total (average) mole fraction of the 1st 
component in the whole surface region, X44 is an analogous 
local value reffered to the homogeneous surface patch with the 
difference of adsorption energies of both components £ 12 7 
£4- Eos but 4 denotes the interval of possible changes in £4 2 
and usually A =(- , + œ) [5], A (E42) is the differential 
distribution function of adsorption sites with regard to the 
value CI- ; 


The distribution functions apperaring in eq. (1) fultils 
the normalization condition: 


J A (E12) 4645 * 1 (2) 


r 
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Equation (1) may be solved analytically by means of the Stieltjes 
transform method [5] or by using the condensation approxination 
approach [s, 25_] „ As the kernel of integral equation in question, 
i.e., as the local adsorption isotherm, the famous Everett’s 
expression is introduced and then eq. (1) may be rewritten in the 
following form: 


+ oo 
xt (40 = fx Xia (K N,) 7A (642)8 8 (3) 


where the constant K,, is defined as follows 14 : 


Kio K°(T) exp (£42/RT) (4) 


Above, k9(T) denotes the temperature dependent constant; it is 
common practice in the theory of adsorption from solutions to 
admit that K (7) 2 1 [14]. 
The general form of the local variable X42,1 corresponds 
to a model supposed for heterogeneous surface as well as to a 
model supposed for adsorption system. The current state of theory 
of liquid adsorption on solids allows us to distinguish four 
models of adsorption systems: 
1. NBP-na model; non-ideal behaviour in both bulk and 
adsorbed phases, where the adsorbed phase is nonautonomous, 
i.e., this phase interacts with bulk one, 
2. NBP-a model with an autonomous surface phase, 
3. IAP model; non-ideal bulk phase and ideal adsorbed phase, 
4. IBP model; both phases ideal. 

Recently, Everett and Podo11 [26] have been 
introduced a model with ideal bulk phase and non-ideal adsorbed 
phase; the suitable abbreviations dealing with this model are 
following: NAF-na and NAP-a. 


168 i *4 . Andrzej Dabrowski 


For IAP and IBP models the variable X421 is independent on 
a model of heterogeneous surface; it means that in this case the 
type of topography of adsorption sites is meaningless [14]. On 
the contrary, both NBP and NAP models require and additional 


e" 


assumption towards the type of heterogeneous surface which 
influences the behaviour of admolecules showing differentiated 
mutual interactions. Hence, for NBP and NAP models also the 
patchwise surface, the variable X, 2,1 has the following form: 


E E 49831 
12,167,140 845 „ for NBP-na model 


Ss: ©) 1 
12, 167,1) 312 for NBP-a model 


“Mat s s AJ EL 1 
21,1% ,1:*1) X42 , for NAP-na model, where = 


T 
or at least 5, 21 


1 À 1 
¥ 21,10% ,11) X45 for NAP-a model, where = 
or at least = 


_ where Ya 2,1/ Na and 83,1 (i=1,2) is the local surface 
activity coefficient; E35 = n/ > ańd p i=1,2 denotes the 
bulk activity coefficients, x] exl/xl and al -/ ag » where 
atex; i (1-1,2) is the activity of the i-th component in the 
bulk phase. 

It follows from the foregoing considerations that for 
patchwise surface equation (1) cannot be solved analytically 
with respect to XP However, if fe more realistic random dis- 
tribution of adsorption sites is introduced, the surface activ- 
ity coefficient fi (i=1,2) appear to be function of the compo- 
sition of the entire adsorbed region and, with respect to model 
of adsorption system used, eventually the composition of the 
bulk phase. Consequently, the variable 1121 corresponding to 


— 
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the random surface is independent on £42 and may be expressed 


as follows: 


for NBP-na model 


s s 1 4 
L21( KU, KK 312 
L246, t) el, » for NBP-a model 


216, b. X1 for NAP-na model 


Wyte Ae RU 1 (6) 
E2058 ,2 X12 for NAP-a model 


1 


* 


822 " for IAP model 
rA » for IBP model 


To take advantage of above-mentioned models with non-ideal 
adsorbed phase, an assumption concerning a solution structure 
is frequently introduced. In this and the simple, statistical 
model of regular solution is applied usually [5,11,15,26]. 
Following by Everett [15], the activity coefficients 
in the regular bulk and adsorbed phases are determined as 


follows: : 
y} = ep -A (1472 (7) 
Mi = exp da (14, t)? + 5 (42 (8) 


where dd /R, Ą”=q”/RT; the parameters a? and q? are respon- 
sible for molecular interactions in the bulk and surface, phase, 
respectively, but at, dss ( , 27. The constants of regular 
solution: 1 and m have their usual meaning and 1+2m=1 [15] 
Equation (8) is valid by assuming that interaction between mole- 
cules in the adsorbed phase equals that in the bulk phase. Tak- 
ing eqs. (7), (8) ana (6) into considerations, ve obtain the 
following expressions for the variable X45: 
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un esp[145( 2x3 10 mą" (2x1-1) ].exp 4"(1-2x7) al, for 
NBP-na model on random surfaces 
exp[14* (2x7 +-1)]exp[A' (-2x) hd], , for NBP-a model 
on random surface 
exp| 14° (2x7 +71 ) + ną" (2x1-1) ]x1, » for NAP-na model 


on random surfaces 


X49" (9) 


exp| 18°( 26 +1) 13, » for NAP-a model on random 
surfaces 

exp[ a? (1-2) kd, + for IAP model 

3 


X45 + for IBP model. 


. Equation (9) will be handled by us during model studies of liquid 


_ adsorption onto heterogeneous solid surfaces. 


ŚĆ THE LF EQUATION AND ITS APPLICATION TO THEORETICAL STUDIES OF 
ADSORPTION PROCESS AT 'SOLID/LIQUID INTERFACE 


To solve analytically eq. (2) with respect to x t* tak- 
, 
ing the local functions (6) or (9) into account, various distri- 
butions x E] 2 ) can be admitted [14]. However, the most popu- 
lar end effective distribution appears to be quasi-gaussian one. 
This function has the following form [3]: 


sin (c) exp (c 8420 


x (£42) = (XRT) (10) 


17 2c Ne) exp, (E 5 exp( A 220) 


where : 

E10 = (£g 7 £42) /RT (11) 
The value of & corresponds to such value of 842 for which the 
. function (10) reaches a maximum. The heterogeneity parameter c 


determines the shape of the distribution (10) and ce(0,1>. In 


i - 
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the case of homogeneous surface, the parameter c=1 and the dis- 


tribution (10) becomes the d -Dirac function for which we have: 
oo 


* * 
J A Cae e je (K- 812) = 1 (12) 
0 co 


Combining eqs. (3) and (10) one gets, by means of Stieltjes 
transform method, the above-mentioned Langmuir-Freundlich 


equation: 


x09) (K x / [14,264 J (13) 
where the variable X42 is given by equations (5) or (9) and the 
constant f, is defined as follows [4]: 


Rio * K(T)exp (E /RT) (14) 


However, the parameter c is temperature dependent [27]: 

be, Seg e=c(T)= c RT (15) 

where the temperature-independent constant c is characteristic 

for a given adsorbent. On the basis of eqs. (9) and (13) we can 

obtain the following expression for. NBP-na adsorption model on 

random surfaces: 

c ln K Z = [65 /r ,)/(xq/e2) Jk ct. (1-n) (2x1-1) + 

(16) 
614 (1-2 +) 

where inf, 2 C/ RT. 

An expression analogous to eq. (16) but valid for homoge- 
neous surfaces was derived previously by Everett [15]. 
Equation (16) enables us to evaluate the isosteric heat of 


adsorption d by using the following formula [28]: 


as (o r) = RT? (Dime J0T) s =(RT*/xpez Yoxz fam) , e (17) 
| 1,t 


" 1,t 


” A y = ^ a la , - 
M es / 
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Bearing in mind eqs. (16) and(17), also applying rules of differ- 
-entiation of the complex functions, we can get the expression 

t 
Qs 


: for » Which corresponds to the most general NBP-na model of 


adsorption system: 
561, 73, 1, 0) = (RTZ xa (l (18) 
where 


| Fgeexp(- £4)(x55 4) “Le, -(1/c) 1n *12, )* xt. X exp(X) (19) 


fà = Lady (g -C) (20) 
x [a (122) (1-23) - 14°(1-2x5 7 ] (21) 
E4 = €,/RT | (22) 


oben for ast, corresponding to all other adsorption 
models, can be developed from eqn (16) and (18)-(22). As an 
instance, for the simplest IBP model of adsorption system, eqs. 
- (16), (18)-(21) become: 


st s 
Q?" = 80 — (1/e,) 1n X12,t 


Equations (16) ana (18)-(22) give the possibility to study nume- 
rically the combined effects of the surface heterogeneity, inter- 
molecular interactions in both surface and bulk phases as well 

as mutual interactions between species belonging to these phases 
on the shape of the excess adsorption istherms and the isosteric 
heats of adsorption. The theoretical adsorption isotherms can. 
de calculated by means of the following equation: 


e 1 8/8 1 
nf a6) = ne (x4 x - x4) (23) 
where nt is the surface phase capacity on the whole heterogeneous 


surface. All model calculations presented below were carried out 


for n£e. 
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Besides, eq. (16) gives a possibility of obtaining criteria i 


for the formation of S-shaped excess isotherms, i.e., an interval 1 


Ld 


of values of Ryo over which an adsorption azeotrope may be observ-. 
ed can be specified, Bearing in mind that at the azeotrope point 

X t - xi, and supposing NBP-na model with ald, m=0,25, 

1=0,5, the following relationship from equation (16) may be 
developed: 


E. - od = exp cdm 2xl-i (24) 


where Rio = (R . Hence, for a given value of 4, an adsorption 
azeotrope is observed in the interval of Rt > given by: 
Kfz E «o, ©) (25) 
Let us compare this result with that arising from homogeneous 
approach presented by Everett [15]. For c: eq. (24) 
becomes to: 
Kyo = exp[ na (2x1 - 1) ] (26 
from which Everett^s condition results [15]: 
Kyo € <exp (-nQ), exp(+m§)7 (27) 
Hence, for the most tightly filled spatial lattice (180,5, 
m=O , 25) a surface azeotrope is observed over the following, 
limited range of values of K45: 
Kjo € «0.606, 1.6487 (28) 
Nevertheless, as far as an adsorption on heterogeneous surfaces 
is considered, the relation (25) is in force. 
Relations analogous to (24), but valid for IAP and IBP 
models of adsorption systems are following: 
ft = ews exp câ} (2x1-4)( 221) ], for IAP model (29) 


and : 
2 = (44 he), for IBP model (30) 


08 09 03 os os 

1. Excess isotneras for adsorption zroa regular binary 
tions Galculated accorcing to eqs. (16)ana (25) by using: 

2 (A) and c=1 (B), §l=05=1, m=0.25, 120.5. Curves (a) - (d) 


note the situation when inf,, runs the values: 5, 2, 0.6, 0, t 


Concluding above resul*s we can state that, in contrast to 


rption on homogeneous surface, the so-called S-shaped iso- 


> value of K,5, which varies from K,,=148,4 (curve a) to 
1 (curve d). On the contrary, in diagram B (c=1, 8 
i E all curves are U-shaped except curve (a) for which 


TS Jae. Here, an adsorption azeotrope is observed as a consequence 


r (4 7*9 hal 


/ 
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of non-ideality of the both phases. What is more, the interaction 
between molecules in the adsorbed and bulk phases is taken here 
into account and the neglection of that is revealed for 1n, 2-2. 
In both parts of Fig. 1 curves (b) denote the situation when 
m=0,25 ( NBP-na model) but dashed lines correspond to m=O ( NBP-a 
model). It follows from Fig 1A that the influence of m is neglé- 
gible during adsorption onto heterogeneous surfaces. On the 

other hand, this effect announces a bigger ascendency when homo- 


geneity of solid is admitted. 


03 0.6 09. 03 0.6 


x żę * 
Fig. 2. Excess 8 for basseption from regular M solu- 
tions calculated according to eqs. (16) and(23) by using: 1nK, 520.6, 
m=0.25, 180.5. Other adsorption parameters are following: (A)c=#0.3, 
and dlz, 45=1 (a); -=- (b); qle-2, 45-1 (c); (B) c=1 and 
gl=2, 89-1 (a); A'=4"=0 (b); &l--2, 49-1 (c); (c) 4-4 and 
c=0.2 (a); c=0.4 (b), c»0.6 (c). In all parts of the figure the 3 
IAS model behaviour (Ąl=4$=0, c=1) is denoted by dotted lines 


The next Fig. 2 shows in parts A and B an influence of the — 


interaction perameters in the bulk and surface phases on the 
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course of the excess adsorption isotherms assuming heterogeneity 
(2A, c=0.3) and homogeneity (2B, c=1) of the adsorbent surface, 
Moreover, in the part C of Fig. 2, the effect of surface hetero- 
geneity on.the excess adsorption isotherm coresponding to IBP 
model is shown. All presented here calculations were carried 
out for 1nK, 2-0. 6, 1. e., for K. 21.82. For comparative purposes 
in all parts of Fig. 2, the dotted lines represent results corres- 
| ponding to so-called IAS model ( ideal adsorption system) for 
which $1-43-0 and c=1. Figure 2 constitutes confirmations of the 
earlier suggestions towards the dominating role of surface hetero- 
geneity in the relation to non-ideality of the liquid solution 
In the controlling the adsorption process in real solid/liquid 
systems. This dominating role of surface heterogeneity will be 
explained in more detail in the next part of present paper. 
It is worthwile to stress here the similarity between 
curves presented in parts A and C of Fig. 2. However, the 
S-shaped isotherms displayed in Fig. 2A vary as a result of 
interactions between molecules of liquid mixture, but these 
showed in Fig. 2C vary as a consequence of changes in surface 
heterogeneity. This fact points out that sometimes an influence 
of surface heterogeneity on the adsorption process at solid/liquid 
interface may be similar to that of the interactions in the liquid 
mixture. 
We will return to this problem in the section dealing with . 
the problem of surface activity coefficients. In the next Fig. 3 


the isosteric heats of adsorption, QSË 


vs xŠ » calculated 

1,t 
according to eqs. (18) - (21) are presented. Part A displays 
the influence of an alternation of the heteregeneity parameter 
c on the heats in question. The curves (a) - (a) correspond to 


the alternation in c from c=0.4 (curve a) to c=1 (curve d), but 


Effects of Surtace Heterogencity in Adsorption os. 177 


03 , 06 09 03 . 06 09 
x M X 1 

Fig. 3. Isosteric heats of adsorption for adsorption from regular 
binary solutions calculated according to eqs. (18) - (22) by 
using: m=0.25, 1=0.5;(A) lnf, 22, d-, 45=0.5 and curves (a) - (a) 
are for: c=0.4 (a) , co. 6 (b), c«0.8 (c), c=1 (a); (B) c=0.7 and 
curves (a) - (d) are for: 41=1.5, $921 and Ink, 4 (a) and 
Ink, 2-2 (c), 4720.5, 4520.5 and Inf, -U (b) and irf, «2(a). In 


both parts of the figures the IAS model behaviour (Ą41=4$=0, 


e=1) for K 2. 2 is denoted by dotted lines 


other adsorption parameters were supposed as follows: Ink, 5=2 
(Ky 2-7. 58), dt, 4520.5, 120.5 and m=0.25. The dotted line 
represents the IAS behaviour ( 4l-49-0, c=1), It appears from 
Fig. JA that the isosteric heats of adsorption from solutions 
are decreasing functions of xi t? however these functions alter 
more rapidly in the case of high heterogeneity (curve a) than 
for homogeneous surface of solid (curve: à). It can be seen, that 
curve d is approximately parallel to dotted line reflecting 
IAS behaviour. It can point out that, in contrast to influence 
of the surface heterogeneity, the shape of functions ast vs „ „t 
is slightly dependent on the values of interaction parameters 
at and e. Such deduction is strongly confirmed by the shape and 
s 


course of the theoretical functions ast vs *, t presented in 


Fig. 3B. All functions showed here were evaluated for c=0.7. 


— 


"d T" 
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Both curves (a) and (b) correspond to NBP-na model of adsorption 


system with 1nf,,=4, however in the former case 41=1.5, &9-1 

(curve a), but in the last case $.-88-0.5 (curve b). Analogous 
calculations, but carried out for lnf,„=2 are represented by 

curves (c) and (a). As in Fig. 3A the dotted line reflects the 

IAS behaviour vith Ink, 22. Obviously, the influence of interac- 
tions parameters on the isosteric heats of adsorption is rather 
greater if the less value of K45 is assumed. Identical conclusion 
towards the influence of interactions parameters on the isosteric 
heats of adsorption from solutions were obtained on the experimental 


grounds by Kiselev and Pavlova [29]. 


SURFACE ACTIVITY COEFFICIENTS RELATING TO THE (1F) EQUATION 


Let us now consider the problem of surface activity coeffi- 


cients in the light of Langmuir-Freundlich equation. First, we 


shall attempt to take advantage of the full form of eq. (16) for 


description of adsorption of benzene from.cyclohexane on fourth 


various adsorbents (see Table 1). The liquid mixture in question 


may be treated as a near-regular solution [50 J and the bulk 


interaction parameter is equal to: $1-0.5 [51]. Let us suppose 
that the parameters â? and 4° appearing in eq. (16) are propor- 
tional to each other. Then 

45/4) = n. 8. (31) - 
Combining eqs. (16) and (31) we obtain the following expression 
(120.5, m=0.25): 

s 1 \-c 1 s 1 1) | [-1 
xi en[1 156) exp[ d 25en4 (1-2x7 .) d 0.75 G- I (32) 
where, according to (24): 
s e s 
X en (f, p/m) + x (55). 
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Equation (32) is the four-parameter one with AA following 
„unknowns: n$, Kiar n and c. To reduce their number, the para- 
Toeter nę: corresponding to system injquestion has been taken 
from the paper by Dabrowski and. --g-ar'o;n e 
[4] where systems from Table 1 vere e bated by means of 
TAP model corresponding to (ir) equation (12). The results of 
bur analysis are given in Table 1. It appears from this Table 
Duas, except adsorption on charcoal, the parameter n is less 
| from unity, i.e., the interactions between admolecules are less 
from those in the bulk phase. Such result is in agreement with 
earlier studies presented by A 8 h et al. [730 and coincides 
with the suggestions of several authors [14]. However, for , 
adsorption on charcoal the molecular interactions in the surface 
Ed appear to be greater than those in the bulk phase (n=1.27), 
but the parameter c has the least of all value among values of 
She characterizing other adsorption systems placed in Table 1. | 
Concluding above results we can state that heterogeneity of the 
solid may either decrease or sometimes increase a non-ideality 
of the adsorbed region. 

_ Equation (32) has not the simple form and its application 
seems to be little encouraging. Besides, this equation does not 
give a possibility to study an influence of both surface hetero- 
-geneity and interactions in the adsorbed phase on the non-ideality 
of the adsorbed solution. This question is very important in . 
liquid adsorption at all because, as it follows from foregoing 
. model studies, either surface heterogeneity or non-ideality of 
“the adsorbed solution may lead to S-Shaped excess adsorption 
| isotherms. What is more, some adsorption systems may be quite 


^ formally interpreted both in terms of the non-ideality of the 
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liquid mixture as well es in terms of heterogeneity of the ad- 
sorbent surface (see Fig. 2). Because the bulk activity coeffi- 
cients are independent on the adsorption experiment it would 
be desirable to give a way of separating of a non-ideal beha- 
viour in the surface region arising from intermolecular inter- 
actions between adsorbed molecules from that due to the surface 


' heterogeneity of solid. 


ne 


0.1 


— — md 


| 

| 
as as as 03 05 09 
Fig. 4. Excess adsorption isotherms (A) and corresponding to 
them functions (34) evaluated for NAP-na model with c=1 (solid 
lines) and for IBP model with c#1 (dashed lines). Other adsorption 
parameters are equal to: Ry 572.718 83-0, 4'=-2, c=1 (solid line 
a); &l-85-0.5, c=1 (solid line b); 47=45=0, c-0.9 (dashed line 
b). In both parts of the figure the IAS model behaviour (AAo, 
c=1) is denoted by dotted lines (c). 


At first, we shall discuss the results of the model studies 
presented in Fig. 4, In part A the solid lines present the 
excess adsorption isotherms calculated according to eqs. (16) 


ek 9, 89-0 (dashed line b). The dotted line (c) presents 


e excess ndsorption BOSE reffering to IAS een (c=1, "EN 


472 2/4, t= rG ), lines (a) and (b) (34) 


* 3 n vua - t(x), dotted line (e) (35) 


n A, IAS denotes thi shonse eee isotherms calculated 


(16) in terms of various, homogeneous models of adsorption 


tems, assuming: 7,,=2.718, n$24, cal, m=0,25, 120.5, and 
b, 12 t 
Vm . 


B 2. these isotherms vere approximated by the following linear 
orm of eq. (12), corresponding to the heterogeneous IBP model of 
| adsorption system (@+=45=0, ): 


Mi de LG t) /x1 Poh =c an (-/- c inf, (56) 


+ xe / f : 
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Table 2. Adsorption parameters of eq. (13) reflecting the behav- 


iour of various models of adsorption systems 


Adsorption parameters of eq. 
system with: 13 evaluated for IBP model 
K45=2-718,c=1, of adsorption system 
1=0.5,m=0.25 


Hence, the parameters c and Rio were obtained and ones are 
given in Table 2. Results presented in Table 2 confirm strongly 
the conclusions obtained on the basis of diagrams presented in 
Fig. 4. Namely, there is a formal possibility to interpret the 
deviations of any adsorption system from IAS behaviour in terms 
of: 

1. non-ideality of one or both phases of adsorption system . 
under condition that cz1 (homogeneous surface of solid), . 

2. effects of the surface heterogeneity of the solid, under 
condition that ĝl=ĝ5=0 (IBP behaviour). 


1 — iF Aubin vd. 


v oe o 4 > - 
- ^ + P — 
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. Obviously, during interpretation most of adsorption systems of 
practical interest both these factors should be taken into 
account. Moreover, it appears from Table 2 that supposition 
about unrealistic values of parameters âl and 4” leads to 
unphysical values of parameters c, which appear to be greater 
| from unity. One must stress, that in the light of the results 
presented in Table 2 and those showed in Fig. 4, the equ poise 
of the NAP-na model with c=1 and IBP model with c#1 is w thout 
controversy. Similar conclusion has been obtained earlie in 
terms of two patch heterogeneous model of the surface by Ou 
and Myers [34]. The above results will be desirable for 
discussion literature adsorption ante interpreted previously by 
Paine a little realistic NAP-na model of adsorption system ( see 
- following). | | 
"RT Let us pass now to a more general considerations dealing 
"with the possibility of separation a non-ideality of the surface 
region caused by intermolecular forces from that generated by 
surface heterogeneity of the solid. In order to do it we shall 
take advantage of results of the osmotic theory of adsorption 
of gas mixtures onto solid surfaces [55] Followed by this 
theory the surface activity coefficient x i (i=1,2) characteriz- 


ing a non-ideality of the adsorbed solution is generated by two 
main factors: 

1. lateral interactions in the surface region, and 

2. heterogeneity of the adsorbent surface. 
For purposes of adsorption from solutions, the former factor 
should be extended over all interactions in the adsorption space, 
i.e., over lateral interactions between adsorbed molecules and 
interactions between molecules belonging to the adjacent bulk 


and surface phases, tóo. 
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Let us suppose now, as previously, the random distribution 
of adsorption sites onto solid. Then, the surface activity coeffi- 
cient 11 may be represented as follows [35]: 


Yi i, n M, int for i=1,2 (37) 
where Vin is the factor responsible for non-ideality of the 
"whole adsorption region, which is generated by the heterogeneity 
of the adsorbent surface, but Ert denotes the factor respon- 
sible for non-ideality of adsorption space caused by the inter- 
molecular interactions. The osmotic theory of adsorption gives 


the following expression for 4. [55]: 


Yin 64, Ue for x5 tÉ (0,1) and i=1,2 (38) 


Keeping in mind the above considerations, the variable * 1 2 
corresponding to homogeneous adsorbent surface has, for the 


NBP=na model, the following form (see eq. 6): 


1 s - 
X42 = 872 + J21,int (59) 
However, an analogous variable, but refering to adsorption on 


“heterogeneous surfaces may be expressed as follows: 


X12,h = 912 * K | (40) 

where, according to eqs. (37) and (38 ) : 
fm = Youn * L21, int (41) 

but Pp 
i. ; L 4, 7$] [(1/e)-1] for xt «€ (0,1) (42) 


Obviously, for c=1 ( homogeneous surface ) X4 2,n*4 2* 
Let us suppose, that adsorption process onto heterogeneous 


surfacesis described by the Everett-type equation [4 5]: 


xi e = Kq2842 J2q/ [442512 12. (43) 
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This equation has the identical methematical form as equation 
relating to adsorption on homogeneous sürféces, However, eq. 

(43) combined with eqs. (41) and (42) leads to the Langmuir- 
-Freundlich isotherms (13), which has here the following form: 


Xi 4 = (f 212 T m) Z LC 204, Us (ts) 


Hence, it is evident that all deviations of any adsorption 
system from IAS behaviour predicted by simples form of Everett’s 
equation: 
x$ K H= K12%12) (45) 
are produced due to: 

1. non-ideality of the bulk phase, 

2. non-ideality of the adsorbed phase, 
but this least factor may be generated by: 

1. intermolecular interactions in the adsorbed region, and 

2. surface heterogeneity of the adsorbent. 

Using eq. (2) we lose a possibility to assess an influence 
of surface heterogeneity on the non-ideality of the adsorbed 
phase and surface activity coefficients corresponding to eq. 
(43) do not give realiable informations about interactions in 
the surface region, because they implicitly include the surface 
heterogeneity effects. 


Obviously, a separation I. between [ 1 p and x is 
U 


s 
i,int 
difficult in general, hovever, such procedure is possible in 
some simple cases (see following). 
At present, we can explain on the basis of foregoing con- 
siderations the similarities arising between NAP model with 
c=1 and IBP model with c#1 (see Fig. 4). For simplicity, the 
NAP-a model of adsorption system will be handled by us for 


which, according to eq. (8) we can get (m=0 ): 


et ; ś . 4 dn UN CN S ur 
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2 
= 0.5 1-xi 46 
li> 25 int exe. T ( P ] for NAP-a Zh Ng, with c=1 ( ) 
(EET „int=exPLO.54 (2x7 ‘on (47) 
However, for IBP model the following relationship is true: 


i - dia LO 
Ka -N, 6 ni 07 (49) 


h 


s rs s 
Table 3. Values of parameters 11 yh? HN, int and Lei, n- Ni, int 
evaluated according to eqs. (46) - (49) for NAP-a model with c=1 


and for IBP model with c=0,9 


NAP=a model with c=1 IBP model with 
and GS=-0.5 c=0.9 


It appears from eq. (43) that the above parameter 154 influence 
the shape of excess adsorption isotherms for NAP and IBP models 
of adsorption systems. 

In the Table 5 the values 12 evaluated according to eqs. 
(47)for 45==0.5 and (49) for c=0.9 and X] „t=0-1> 0425.5 | 


are listed. It can be seen that these values are practically 


E 
or IBP model with r! (48) 


^ 
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Xu xy 
Fig. 5. Surface activity coefficients J vs X] t (A) and corre- 
sponding to them functions r 21 vs x4 t calculated according to 

?, 


eqs. (46)-(49) for O. 7, 81 


equal in the whole range of NT thus in the whole range of 


xi. In view of the small values of excess adsorption isotherms 


for > t (and x) approaching to unity, the differences between 
, 


1 
1 


are observed in Fig. 4. In the Table 3 the values of 11 both 


suitable functions (34) in the region of high values of x 


|. for NAP-a model with c=1, 4”=-0.5 and IBP model with c=0.9 are 
_ listed, too. It follows from Table 3 that the contribution to 


the coefficient H- Rn . Hi, int from In is like, for 
a? <0, to that arising from n any * One of possible consequences 
, 


of this fact is presented in Fig. 5, where the model studies- 
carried out on the basis of eqs. (37), (41) ana (46) - (49) for 
c=0.7 and $9241 (positive deviations from ideality) are display- 

ed. Part A in Fig, 5 presents the following functions: Yin vs Xi t 


s s 8 8 
Ji „int VS X4 at and {i= Yin H, int vs X4 t · Irrespective of 
the fact that 5-1, the surface activity coefficient 1 HAST Y 


‘phy vi. 


from unity in the whole concentration interval. It means, that 
due to surface heterogeneity, an adsorbed region is characteriz- 
ed by the surface activity coefficients Ü i suitable for liquid * 
mixture showing negative deviations from ideality. In part B of 
Fig. 5,the functions T2 vs X] t’ (> jh Ys Xi t and Li, int : 
vs xit are shown. It follows from diagram B that over concen- 
tration interval x4 tÉ (0.0.5) the influence of 15; jh on the 
shape of excess adsorption isotherms is greater than that. 
arising from the factor 421, int however, for x4 E (0.5, 1) 
the situation is: inversed. Such conclusion is strongly supported 
by the results given in Table 5. 

Let us now consider the utilization of eqs. (37), (41) and 
(46) - (49) to characterize two literature experimental systems. 
In order to do it the parameters nl, c and 8? have to be obtain- 
ed on the basis of (IF) equation. In the paper [5 Da brow- 
s k 1 and Jaroniec have been applied eq. (13) in 
the form corresponding to NBP-ra model of adsorption system, for 
which the bulk phase was admitted to be non-ideal (but not 
regular ) and the surface phase was supposed to be regular one 
with the autonomous surface phase. Supposition about the autos 
omous surface phase, Supposition about the autonomous character 
of surface phase seems to be reasonable in the light of model 
studies presented in Fig. 1. 

First, we shall discuss the results of our analysis carried 
out for the system: piridine (1)+ethyl alcohol (2) on charcoal 
at 293K investigated experimentally by the B lackburb 
et al. [32]. Above system is characterized by the following set 
of adsorption. parameters [5]: n£25.58nnol/g. c=0.84, G”=0.915 
and a0 [32]. The surface mole fraction appearing in eqs. (46)- 


— è 


03 06 09 i 03 06 09 


piridine (1) + ethyl alcohol (2) on charcoal at 293K [32]. con- 
tributions to the surface activity coefficient I 1 (solia line 
with circles) (a) and to corresponding to it parameter + 24 
(solid line with circles) (5) from the interactions in ihe 
adsorbed phase (strongiy dashed curves ) and surface hetero- 
geneity effects (dashed curves ). In the diagram (4) the bulk 
activity coefficients are represented.by solid line without 
les. In both parts of the figure the dotted lines calculat- 
ed according to the equations representing I Ti ERT r dat 
and Jai, [24,n* J Zl, int by using: c=0.84 and 4*=0.915 [5] 


-(49) has been obtained according to eq. (33). to avoid the 
1appi f s s rs s 
overlapping of the functions L 1 vs X ot and F 2 VS X4 yt? in 


Fig. 6A the former one is presented only. The solid line parallel 


to concentration axis, with ordinate equal to unity, presents 
the bulk activity coefficients: B =1 (i=1,2), but the solid. 
line with circles experimental points shows the function 


{i= r 1 (A, t) · It can be seen that the experimental results 


€ region aud the otted lines n 
Me present the courses of suitable functions evaluated by means 
of eqs. (46)-(49) for lower values of x1 ,t* In spite of the 
FE that the solution is ideal in the bulk phase, the functio: 


x i tha ae 
not cover the ow" x 


m ffs i, int VS x at points out an existence of lateral interactions R 
between admolecules caused probably by the mainly chemical  - " 


unhomogeneity of the charcoal surface. However, due to course oa 


LI 


7" eoefficient:: ys. 42 15 „h * 15 int 18 approximately equals to ZY Y „(3 
WSZ H 


" s 
ir 
In the part B of Fig. 6, the function 15 vs CM z ia 
A I 

"shown, too. In the next Fig. 7 the results of analogous dn. 15 | 
* 


NA 228 for the system: n-heptane (1 )+cyclohexane (2) on tin-oxide ge 


~ 


at 308K are presented. This system was investigated exper iia 


of adsorption a ERINA ek nę=0.69mmo1/g, c=0,93 and 85-0. a 
In Fig. 7A the both functions ¥ 4 (11, 2) vs Xi t are showed 
The solid line without circles presents the bulk activity coe- 


the, function ri 1 vs x, t goes below the function J 1% A 
and is greater from unity, In Fig. 7B the function rs 21 Vs 11, 
is presented and this one has different course from that showed 
in Fig. 6B. It follows from Figs. 6B and 7B that the dominating — 


role of the surface heterogeneity on the shape of excess AA 


tion isotherms deals with the lower values of x, t but in the 
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Fig. 7. Surface. activity coefficients for the experimental 
system: n-heptane (1 )+cyclohexane (2) on tin-oxide gel at 308K 
[56]. Contributions to the surface activity coefficients fs 
(solid line with white circles) and {3 (solid line with black 
circles) (a) and to corresponding parameter In (solid line 
with white circles) (B) from the interactions in the adsorbed 
phase strongly dashed curves and surface heterogeneity 
effects (dashed curves). In the diagram A the bulk activity 
coefficient 5 is represented by the solid line without circles. 
In both parts of the figure the dotted lines calculated accord= 
ing to the equations representing I (in T int (i=1,2) and 
I2 Li, n= J21,int by using: c=0.93, 4° = 0.838 [5] 


region of higher mole fractions, the dominating role of lateral 
interactions between admolecules is observed. 

Finally, we can explain a reason for which the S-shaped 
isotherms are produced by (17) equation on strongly heterogeneous 
surfaces, irrespective of values addmitted for Rio and other 
adsorption parameters (see Figs. 1 and 2). In Table 4 the 
values of function Y 2n VS C, obtained according to eq. (42)- 


| 
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Table 4. Values of VERS vs c evaluated according to eq. (49) 
s s 
for x4 {70.1 and for x4 ,t70+95 


Fian = Lein (9) for 


0,72 
0.48 
0.28 
0.14 
0.05 
0.01 
14.04.1079 
6555.11 7.80.1075 


8 


3.87.10 3.23.1012 


for xà „t=0+1 and x3 £70-95 are given. It appears from Table / 

that the difference betveen In ,h at 4. and x4 420695 is 
insignificant, but for c=0.1 this difference is immense. Hence, 
bearing in mind eqs. (33) and (43) we can state that the excess 
isotherm is greater from zero for low values of X] t (thus x1) 
and less from zero for higher values of xi t (thus xi). Due to 
this fact an adsorption azeotrope appears on strongly hetero- 


geneous surface of solid. 


— > 


DATA MEASURED AT VARIOUS TEMPERATURES 


de | x 
5 Recently Everett and Podoll _26] have 


b 


ee 


tal excess adsorption isotherms circles in comparison to 

e theoretical curves solid lines calculated according to eq. 
51. In all parts of the figure the dashed lines present the 

io „7 RAN isotherms calculated according to eq. 51 with 

A by means of the parameters L. and nę from Table 5 


© 


alysed the excess adsorption isotherms measured for adsorption 


ie 8. Adsorption of benzene 1 from ANTE eet m 2 on 
G caphon at 273K A, 298K B, 313K © and 333K D 26. Exper- 


15 3 : | 
infh xL Me] i^ 


Fig. 9. Adsorption of benzene 1 from 1, „Z=dichioroetianę 2 
on Graphon at LO temperatures 26 , linear dependence e 
| calculated PARE to eq. 26. Circles denote experimental | 


points 


196 ^ Andrzej Dabrowski 


positive deviations from Raoult’s law Be]). On the other hand, 

! if follows from above Presented considerations that no possibility 
exists to interpret S-shaped isotherms in terms of IBP model of 
adsorption system with c=1 ( homogeneous surface). For this reason, 
to explain experimental results by means of homogeneous approach, 
Everett and Podoll [26 ] aanittea the NAP-na model 
of adsorption system with large negative deviations from ideality 
in the regular adsorbed phase. Such same model has been utilized 
by. L. i and Ku [39] to explain experimental results dealing 
with the adsorption of some bulk near-ideal liquid mixtures on 
silica gel surface, Such approach was severely criticized by 

^a MZC ar [12] who stated that admission of large deviations 
from ideality in the adsorbed phase is little realistic in view 
of the fact that the bulk solutions are nearly ideal and show 
certain positive deviations from Raoult”s law, On the other hand, 
the formal equipoise between NAP.model with c=1 and IBP model 
with c1 has been proved and explained by us. For this reason 

"we will try to explain experimental results measured by 
Everett and Podoll [26] by means of IBP model 

of adsorption system, assuming H. In order to do it, eq. (36) 
has been utilized and results of our numerical analysis are 

| given in Table 5. During our calculations the values of Xi t 
appearing in eq. (36) were evaluated on the basis of eq. (33) 


where the parameter nę was estimated as follows: 
nę = a / ag = 0.365 mmol/g (50) 


Above, the specific surface area of the Graphon, as, Was assess- 
ed by the authors of the experimental data and a =84m?/g [2], 


but as denotes an arithmetic average of the cross-sectional 
ur / 
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Table 5. Adsorption parameters of eq. (56) evaluated for adsorp- 
tion of benzene. (1) from 1,2-dichloroethane (2) on Graphon at 


various temperatures 


areas of the components of liquid solution and, following by 
Everett and Podoll, 922.315 10? m?/umo1. 

Figures 8 and 9 show the results of our analysis carried 
out by means of eqs. (33) and (36). Parts A, B, C and D of 
Fig. 8 present experimental excess adsorption isotherms 
(circles) measured at 273K, 298K, 313K and 323K, respectively 
with comparison to the theoretical curves solid lines calculat- 
ed by assuming the best-fit parameters from Table 5, by means 
of the following formula Ca]: - 

nf -n ada (4-xt) cf] / Kiez) È (51) 
However, in Fig. 9 the linear e calculated according 
to eqs. (33) and (36) with comparison to experimental points 
(circles) are presented. It appears from Figs. 8 and 9 that the 
experimental data in question, measured at various temperatures, 
can be described relatively well by the () equation correspond- 
ing to IBP model with cZ1. The dashed line in each part of Fig. 8 
shows the corresponding homogeneous surface excess isotherms 


evaluated according to eq. (51) with ©=1 and using the same 


x rs 
N 


7 d i 
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ma 29 M M ry 
Fig. 10, 1R} vs(1/T) (A) (solid line with black circles) and 


c vs T(B) (solid line with white circles) for adsorption of 

benzene (1) from 1,2-dichloroethane (2) on Graphon at various 

temperatures [26]. Values of suitable parameters are given in 
Table 5 


values of 42 from Table 5, In Fig 8A considerable influence of 
inappreciable surface heterogeneity of Graphon (see Table 5) on 
the shape of excess isotherms is ebserved. This fact is connect- 
ed with the small values of K 2 corresponding to systems in 
question. It follows from Table 5 that both heterogeneity para- 
meter c as well as the constant 242 are tempera ture- dependent. 
In part B of Fig. 10 the variation of c with T is shown; the 
relation is linear and from the slope of this line the parameter 
cg can be easily evaluated. On the basis of eq. (15) we can 
write: 
c, = (1/R) (ac/az) (52) 

and it has been found that c,=1.75 107% mol/J. 

In contrast to the parameter c, the parameter co is tempe- 
rature=independent and this one characterizes the surface hetero= 


geneity of the Graphon. 


AJ SE 


ink, = Ink e(r) + /RT (55 


(Hence, the slope of the straignt line in Fig 10A directly deter 
bee parameter Ink°(T) and 1nk?(T)-0.144, ices, K?(T)=1.1 


In view of this fact that K (T)is slightly different from unity f 
"v < 


the interesection point of the dz o line with ordinate sa 


Ve approximately equal to: 
| Cain, / (1)71] = Eve 
here it has been found that Co 0.214 KJ/mol=0.057kcal/mol. - 4 


' Introducing the above evaluated parameter R to eq. (10), the a 


> distribution function characterizing the adsorption systems in 
gu question may be easily evaluated. In view of the temperatura "M 
PI» 4 


e of the parameter c, these functions will have a. 


|. fixed value of C O. 2 KJ/mol. 


Table 6. Adsorption parameters of eq. (36) for experimental | V 
systems investigated at 293K on silica gel by Li and Ku [29 


L iquid* wisxtur-rte 


Toluene (1) + Chlorobenzene (2) 0.856 
0.898 
0.917 
Benzene (1)+ Chlorobenzene (2) 0.973 


1 + Bromobenzene (2) 0.982 


Toluene (1) + Bromobenzene (2) 


Benzene (1) + Toluene (2) 


Chlorobenzene 


a) Obtained by using eq. 50 by means of parameters a. and S 
x 
3 from the paper [ 39). 
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Fig. 11. Linear dependences calculated according to eq. (36) 
for adsorption systems from Table 6 [39:] (a) toluene (1) « chlo- 
robenzene (2), (b) toluene (1) + bromobenzene (2), (c) benzene 

(1) * toluene (2), (à) benzene (1)+ chlorobenzene (2), (e) 
chlorobenzene +  bromobenzene (2). Circles donote experimental 


points.. 


Concluding foregoing considerations we can state that little 
realistic NAP model with 4*<o, used in previous analysis of 
adsorption of benzene from solutions in 1,2-dichloroethane on 

-` Graphon at various temperatures can be explained in terms of 

more realistic IBP model of adsorption on a heterogeneous surface. 
This remark is in force also towards experimental systems inves- 
tigated by Li and R u [L 59]. The results of numerical analysis 
dealing with their data are given in Table 6, but in Fig. 11 the 
linear dependences plotted according to eq. (36) by means of 


parameters from this table.are presented, too. 


E t * * 
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STRESZCZENIE 


Stosujac model monowarstwowej fazy powierzchniowej opisano 
proces adsorpcji z dwuskładnikowych roztworów nieelektrolitów 
na energetycznie niejednorodnych powierzchniach adsorbentów 
stałych, W opisie uwzględniono nieidealność fazy powierzchnio- 
wej i objętościowej oraz oddziaływania pomiędzy cząsteczkami 
należącymi do tych faz, 


Peswme 


Menonbsya uonenb MOHOCJIOA AIA AHCOPOLMCHHOŃ ash, ONUCAHO 
nponecc a4copórMu Ka ÓWHADHHX pacTBODOB He3JIeKTpOJHTOB Ha 
SHepreTMWeCKM DETEDOTEHHHX TBEPANX AHCODÓEHTAX. JTO ONUCAHKE 
YUTHBAST HeMJISSJUIBHOCTb AHCOPÓWDOBAHHOŃ M ÓJNOUHOŃ asu M BZA- 
uMOJeHCTBUe MeXjly MOJeKyJIaMM ƏTHX das. 


